Steering Plasmodium with light: 
Dynamical programming of Physarum 

machine 



Andrew Adamatzky 

University of the West of England Bristol BS16 1QY United Kingdom 
andrew.adamatzky@uwe.ac.uk 

(N 

3 

<^ Abstract 

A plasmodium of Physarum polycephalum is a very large cell visible by unaided 
eye. The plasmodium is capable for distributed sensing, parallel information pro- 
Pi^ cessing, and decentralized optimization. It is an ideal substrate for future and emerg- 
ing bio-computing devices. We study space-time dynamics of plasmodium reactiom 
to localised illumination, and provide analogies between propagating plasmodium 
and travelling wave-fragments in excitable media. We show how plasmodium-based 
computing devices can be precisely controlled and shaped by planar domains of 
I illumination. 
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Q\ 1 Introduction 

O 

At one phase of its life-cycle Physarum polycephalum\^\ lives as a huge sin- 
gle cell with many diploid nuclei. Such a life-form is called plasmodium. The 
plasmodium feeds on bacteria, spores and other microbial creatures. When 
foraging for its food the plasmodium propagates towards sources of food par- 
ticles, surrounds them, secretes enzymes and digests the food. Typically the 
plasmodium forms a congregation of protoplasm in a food source it occupies. 
When several sources of nutrients are scattered in the Plasmodium's range, 
the plasmodium forms a network of protoplasmic tubes connecting the masses 
of protoplasm at the food sources. 

Nakagaki et al [20f21|l22] showed that the topology of the Plasmodium's proto- 
plasmic network optimizes the Plasmodium's harvesting on distributed sources 
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of nutrients and makes more efficient flow and transport of intra-cellular com- 
ponents. The plasmodium is considered as a parallel computing substrate com- 
plementary [3] to existing massive-parallel reaction-diffusion chemical proces- 
sors [2J. The plasmodium functions as a parallel amorphous computer with 
parallel inputs. Such a parallel biological computer takes data in a form of 
spatially distributed sources of nutriets, and represents results of computation 
by configuration of its entire body and protoplasmic veins. The plasmodium is 
capable for a approximation of shortest path, computation of planar proximity 
graphs and plane tessellations, primitive memory and decision-making. 

In j3] we shown that the plasmodium of Physarum polycephalum is a gen- 
eral purpose computing machine. This is because plasmodium implements the 
Kolmogorov-Uspensky (KUM) machine [TTfSE] in its foraging behaviour. The 
KUM is a mathematical machine in which the storage structure is an irreg- 
ular graph. The KUM is a forerunner and direct 'ancestor' of Knuth's link- 
ing automata |16j . Tarjan's reference machine [33], and Schonhage's storage 
modification machines [26127] . The storage modification machines are basic 
architectures for random access machines, which are the basic architecture of 
modern-day computers. The plasmodium-based implementation of KUM [I] 
is a first-ever biological prototype of a general purpose computer. 

The key component of the KUM is an active zone [TTll3~6] , which may be seen as 
a computational equivalent to the head in a Turing machine. Physical control 
of the active zone is of utmost importance because it determines functionality 
of the biological storage modification machine. In present paper we experi- 
mentally demonstrate that propagation of KUM's active zone can be tuned 
by localized domains of illumination. 

The plasmodium of Physarum exhibits negative phototaxis. General under- 
standing of the plasmodium response to light is that the plasmodium moves 
away from light when it can and switches to another phase of its life cycle 
or undergoes phragmentation when it could not escape light. If plasmodium, 
particularly the starving one [T3], is subjected to high intensity of light the 
plasmodium turns into a sporulation phase [25J. There are evidences that 
phytochromes are involed in the light-induced sporulation [32] and sporula- 
tion morphogene is transferred by protoplasmic streams to all parts of the 
Plasmodium [TJ]. 

Photofragmentation is another physiological response to strong and unavoid- 
able illumination. When plasmodium is illuminated by ultraviolet or blue 
monochromatic light, in hostile environment of laboratory conditions, it breaks 
into many equally sized fragments (each fragment contains around eight nu- 
clei) [15J. The fragmentation is transient and after some time the fragments 
merge back into a fully functional plasmodium. 
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In the paper we focus on photomovement, the less drastic response to illumi- 
nation than two scenarios mentioned above. Plasmodium of Physarum poly- 
cephalum exhibits most pronounced negative phototaxis to blue and white 
light [91130] . The illumination increase causes changes in the plasmodium os- 
cillatory activity, degree of changes is proportional to a distance from the 
light source [37ffTTj . Exact mechanisms of the response to light is yet un- 
known. There is however a few phenomena uncovered in experiments. First 
is presence of phytochrome-like pigments [T5] , which might be primary recep- 
tors of illumination. The light-response of the pigments triggers a chain of 
biochemical process [29]. These processes include increase in activity of iso- 
merase enzymes |31j . changes in mytochondrial respiration [18], and spatially 
distributed oscillations in ATP concentrations [35] . 

Nakagaki et al [24.23] undertook first ever experiments on shaping Plasmod- 
ium behavior with illumination. They discovered that protoplasm streaming 
oscillations of plasmodium can be tuned by, or relatively synchronized with, 
periodic illumination [21]. They also demonstrated that plasmodium optimizes 
its protoplasmic network structure in the field with heterogeneous illumina- 
tion [23J: thickness of protoplasmic tubes in illuminated areas are less then 
thickness of tubes in shaded areas [23J. These indicate that illumination gra- 
dients could be a convenient tool to input instructions to different parts of 
Physarum machines in parallel. Several basic questions need to be answered. 
What exactly a pseudopodium or a migrating Physarum do when approach 
an illuminated domain? How light can be used to program plasmodium move- 
ments? What types of plasmodium reflections can be implemented using light- 
mirrors? Can we shape the Plasmodium's network structure by heterogeneous 
illumination? In the paper we present experimental and theoretical findings 
we obtained while trying to answer the questions. 



The plasmodia of Physarum polycephalum^\were cultured on wet paper towels, 

fed with oat flakes, and moistened regularly. We subcultured the plasmodium 
every 5-7 days. 

Experiments on controlling the plasmodium with light were undertaken in 
standard Petri dishes, 9 cm in diameter. A substrate was 2% agar gel. Light- 
obstacles, or illumination domains, were implemented using electro-luminescent 
sheets Based on previous works in negative phototaxis of Physarum we 

2 Culture of Physarum polycephalum was kindly supplied by Dr. Soichiro Tsuda 

3 Manufacturer Seikosha, supplier RS Components Ltd Birchington Road, Corby, 
Northants, NN17 9RS, UK 



2 Methods 




3 




I ■.■II i Ji<h v. ill'. ikuij 45m v4 iiimI 

* tpw|u. | Jik' i\ ilti 1r;iruipnmiil 
-. 1 1 . ■ ; ■ ..- ..hi Hi 

i kcSi !'-Liiiliiiih:s-.(.jiI sJkuI 



(a) 




Fig. 1. Experimental setup: (a) scheme, (b) photographs of illuminating domains, 
geometrical shapes. 

have choosen blue illuminating sheets. The nominal sheets brightness (at 110 
V, 400 Hz AC) is 73 cd/m 2 (blue). The electro- luminiscent sheets did not pro- 
dure heat, therefore can be positioned in close contact with growing substrate. 
Masks were prepared from black plastic: rectangles and triangles were cut in 
the plastic. When a mask is placed on top of the electro-luminiscent sheet, the 
light is passed only throught the cuts (Fig. \U, 



The experiments were conducted in a room with diffusive light of 3-5 cd/m, 
22°C temperature. In each experiment an oat flake colonized by the Plas- 
modium was placed on one side of a Petri dish, and few oat flakes without 
Plasmodium at the opposite side of the Petri dish (to shape directional prop- 
agation of the plasmodium). Petri dishes with plasmodium were periodically 
scanned on a standard HP scanner. The only editing done to scanned images 
is colour enhancement: increase of saturation and contrast. 

When plasmodium is cultivated on a non-nutrient substrate a profile of its 
pseudopodium's propagating tip, and particularly the profile of the plasmod- 
ium propagating as a whole, is isomorphic to shapes of and behaves as wave- 
fragments in sub-excitable reaction-diffusion chemical media, see details in [5|6] . 

Pictures quality is dramatically reduced. See version published in the Journal for 
high-resolution pictures. 
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Therefore we simulated propagation of the plasmodium with two-variable 
Oregonator equation i 12~ll3~Tl. used originally to simulate light-sensitive Belousov- 
Zhabotinsky reaction with applied gradients of illumination [8~|IT9~] : 




In terms of Belousov-Zhabotisnky reaction the variables u and v represent 
local concentrations of bromous acid HBr02 and the oxidized form of the cat- 
alyst ruthenium Ru(III), e sets up a ratio of time scale of variables u and v, q 
is a scaling parameter depending on reaction rates, / is a stoichiometric coef- 
ficient, is a light-induced bromide production rate proportional to intensity 
of illumination. The <ft is an excitability parameter. Moderate intensity of light 
will facilitate excitation process, higher intensity will produce excessive quan- 
tities of bromide which suppresses the reaction. There is no diffusion term for 
v because we assume the catalyst is immobilized. 

To integrate the system we used Euler method with five-node Laplasian oper- 
ator, null boundary conditions, time step At = 5 • 10~ 3 and grid point spacing 
Ax = 0.25, with the following parameters: <fi = <po — a/2, A = 0.0011109, 
0o = 0.0804, e = 0.03, / = 1.4, q = 0.022. 

Illuminated domain (light obstacle) S of the experimental space was simu- 
lated by higher values of the parameter <p: if given point belongs to S then 
cf) = 0.085 (the medium inside the light obstacle becomes non-excitable) oth- 
erwise <p = 0.0804. For such set of parameters the model represents a sub- 
excitable medium (at the edge between of non-excitability and full excitabil- 
ity). The sub-excitable media exhibits self-localized excitations. A local distur- 
bance leads to formation of traveling excitation wave-fragments which preserve 
their shape for finite period of time and then either expand or collapse. The 
traveling wave-fragments imitate propagating plasmodium. 



3 Trees and waves 

There are two distinct forms of the plasmodium: a protoplasmic tree and 
a traveling localizations. The protoplasmic tree is formed when plasmodium 
forages the space by sprouting pseudopodia in various directions. The pseu- 
dopodia remain connected to original location (and still main "body" ) of the 
plasmodium by protoplasmic tubes. In certain situations (exact characteris- 
tics of such situations might be a subject of further studies) the plasmodium 
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(a) t = h (b) t = 3 h (c) t = 6 h 
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(d) t = 11 h 


(e) t = 20 h 


(f) t = 29 h 




(g) t = 33 h (h) t = 39 h 

Fig. 2. Transformation of plasmodium from protoplasmic tree to travelling localiza- 
tion (migrating plasmodium). 

leaves its original location and propagates as a whole on the substrate. When 
the plasmodium propagates as a whole (migrating plasmodium) it has shape 
of a wave-fragment traveling in the a sub-excitable medium. The boundary 
between two morphologies may be fuzzy, and very often we can observe wave- 
fragment like shapes of pseudopodia, and gradual transitions between tree-like 
and wave-fragment like morphology. 

An example of the "tree to wave-fragment" transformation is provided in 
Fig. [2j An oat flake colonized by plasmodium is placed in south part of a Petri 
dish (Fig. [2^i). The plasmodium sprouts several pseudopodia exploring the 
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space (Fig.[2p,). The protoplasmic branches die out when encounter illuminated 
domains (Fig. ^pc). A group of branching pseudopodia tries to find a way 
around the triangular light obstacle (Fig. [2^) but eventually abandons the 
attempt (Fig. gi). 

By that time bacteria on the original oat flakes are exhausted and the Plas- 
modium switches to its migration phase (Fig. [2]i). The Plasmodium abandons 
its original oat flake and starts propagating as a whole along the eastern wall of 
the Petri dish (Fig. |2kf). A typical wave-fragment of the propagating Plasmod- 
ium is formed (FigH2t) which heads toward the source of nutrients (a group 
of oat flakes in the northern part of Petri dish). Eventually the plasmodium 
reaches new source of nutrients (Fig. [2}i). 

When the plasmodium propagates as a whole it looks like and behaves like a 
wave- fragment in sub-excitable medium 06]. Moreover, migrating plasmod- 
ium is more sensitive to potential environment threats than just propagating 
pseudopodia. This is because a protoplasm in pseudopodia can always be re- 
tracted back to main body of plasmodium, thus pseudopodia can "take risks" . 
The migrating plasmodium is more vulnerable, because any "mis-calculation" 
in choosing its migration route may lead to disaster. 



4 Diverting plasmodium 



In an ideal situation the plasmodium propagates as wave-fragment and — un- 
less encounters an obstacle on the way — keeps its shape and velocity vector 
conserved. If a proximal part of the plasmodium wave comes upon a highly illu- 
minated domain, frequency of protoplasm oscillations in this domain increases. 
Due to a difference in the protoplasm oscillation frequency the plasmodium 
wave slightly turns to the side with less oscillating protoplasm. 

Finding 1 One can steer propagating plasmodium using light obstacles. 

A deflection of the plasmodium wave by light triangle is demonstrated in 
Fig. |3} The plasmodium wave propag ates North-North-East (Fig. |jie). The 
plasmodium hits a light triangle with its western side (Fig. |3]Df). The light 
increases frequency of oscillations in the illuminated part of the plasmodium. 
The plasmodium wave turns eastward (Fig. and travels in the new direc- 
tion until hits the dish's wall (Fig. ^dh). 

The plasmodium may not reflect from light obstacle immediately. After en- 
countering even highly illuminated area the plasmodium wave carries on trav- 
eling in the original direction and only after passing obstacle starts diverting. 
Such example of delayed reaction to illumination is shown in Fig. |4j Plas- 
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(a) t = h (b) t =1 h (c) t =5 h 




(d) t =11 h (e) t =0 h (f) t =1 h 




(g) t =5 h (h) t =11 h 

Fig. 3. Deflection of a plasmodium wave by light triangle, (a)-(d) photographs of 
the plasmodium, (e)-(h) schemes of the propagating plasmodium waves. 

modium wave, originated from the plasmodium colonizing a group of flakes in 
northern part of the dish, travels South (Fig. [4£ic). The wave passes through 
rectangular illuminated area without any immediate reaction (Fig.|4^i). The re- 
sponse occurs several hours later. The plasmodium steers South- West (Fig.^d 
Note that the plasmodium is not diverted by second (lying southward) light 
rectangle, because when diversion happens the plasmodium is far from the 
second light rectangle. 

Another example of the steering by light obstacles is shown in Fig. |5j The 
plasmodium wave, heading North, enters the illuminated rectangular area 
(Fig. |5)ic). The wave is slightly displayed towards eastern side of the light 
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(a) t = h (b) t = 7 h 




(c) t = h (d) t = 7 h 

Fig. 4. Deflection of plasmodium wave by light rectangle, (a)-(d) photographs of 
the plasmodium, (e)-(h) schemes of the propagating plasmodium waves. 

rectangle. Due to differences in light-induced oscillation of protoplasm the 
plasmodium turns North-East and continues traveling in this direction till 
collides with a wall of the Petri dish (Fig. [^d). 

Finding 2 Combining light obstacles and sources of nutrients (chemo-attractants) 
one can implement precise control of plasmodium waves. 



The plasmodium is attracted to oat flakes populated with bacteria: we exper- 
imentally demonstrated that the plasmodium senses and reacts to sources of 
nutrients placed as far as 3-4 cm away of the plasmodium. The plasmodium 
'ascents' along gradients of the attractants till it reaches sources of nutrients. 
By placing light obstacles in the attractive field we can tune and shape trajec- 
tory of the plasmodium motion. For example, the plasmodium shown in Fig. [6] 
is reflected eastward by first illuminated rectangle it encounters (Fig. [6^d). 
Later the plasmodium turns North being attracted to oat flakes in Northern 
part of the Petri dish. The plasmodium is directed westward by the second 
illuminated rectangle (Fig. ^jpe). However it does not continue West but turns 
North again to reach the sources of nutrients (Fig. [6fcf). 
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(a) t = h (b) t = 12 h 




(c) t = h (d) t = 12 h 



Fig. 5. Deflection of plasmodium wave by light rectangle, (a)-(d) photographs of 
the plasmodium, (e)-(h) schemes of the propagating plasmodium waves. 



5 Multiplying plasmodium waves 



Finding 3 A propagating plasmodium wave or a pseudopodium can be split 
by suitably shaped domain of illumination. 

In some situation propagating plasmodium hits a light obstacle which is small 
enough to divert the whole plasmodium wave. If parts of the plasmodium 
wave remain outside the illuminated shape, these parts continue travel as 
"independent" plasmodium waves. Thus the plasmodium wave will be split 
into two waves. An example is provided in Fig. [7| In this particular experiment 
the plasmodium wave is comprised of several pseudopodia moving together in 
a close group. On its way toward source of nutrients the wave runs across 
illuminated triangle (Fig. [7^d). The pseudopodia try to steer away from the 
source of light and move toward western and eastern sides of the illuminated 
triangle (Fig.[7pe). Eventually two separate groups of pseudopodia are formed, 
one group travels North- West, second group moves North-East (Fig.FTbe). 
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(a) t = h 



(b) t = 12 h 



(c) f = 16 h 




(d) t = h (e) t = 12 h (f) £ = 16 h 

Fig. 6. Steering plasmodium with combination of light obstacles and attracting field. 
6 Optimization of foraging in presence of illumination-obstacles 



Finding 4 When Physarum scouts a space with pseudopodia it abandons the 
pseudopodia which encounter illuminated domains and sprouts new pseudopo- 
dia in the less illuminated areas and towards the sources of nutrients. 

Nakagaki et al [23] already demonstrated that the plasmodium optimizes its 
protoplasmic network in presence of light, namely the plasmodium reduces 
number and size of protoplasmic tubes exposed to light. In present paper we 
do not focus on optimization and do not provide statistical evidences but 
rather few morphological examples of how growing protoplasmic network is 
shaped by illuminated domains. 

Most typical example is shown in Fig. [8j Initially an oat flake colonized by 
the plasmodium is placed in a central part of the Petri dish (Fig. [8^c). Three 
intact oat flakes are placed near the Petri dish wall so they form vertexes of 
a triangle centered in the initial position of the plasmodium. The plasmod- 
ium sprouts several pseudopodia heading towards the oat flakes. Most of the 
pseudopodia encounter illuminated areas and seize propagation. The pseu- 
dopodium propagating East-South-East manages to avoid illuminated areas 
and eventually reaches North-East oat flake (Fig. ^pd). 
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(a) h 



(b) 6 h 



(c) 11 h 




(d) Oh (e) 6 h (f) 11 h 

Fig. 7. Splitting plasmodium wave with illuminated triangle. 



There are evidences that the plasmodium may "evaluate" size of the light 
obstacles by sprouting pseudopodia in several directions, integrating informa- 
tion about size of obstacles (possible via protoplasm oscillation frequencies) 
and then making a decision about where to sprout further pseudopodia. An 
example of such decision making is provided in Fig. |9j In the experiment re- 
ported the plasmodium was placed in the center of a Petri dish, and additional 
oat flakes were placed in North and South poles of the dish. The plasmodium 
sprouts several pseudopodia to explore the space. Some of the pseudopodia 
encounter illuminated areas. The illuminated area south of the plasmodium 
is larger than the illuminated area north of the plasmodium. Therefore the 
Plasmodium does not sprout any more pseudopodia southward but produces 
a protoplasmic branch growing North. This pseudopodium reaches the south 
source of nutrients (Fig. [9J. 



Such mechanism of scouting allow plasmodium to pre-shape its foraging net- 
work. When eventually plasmodium spans all sources of nutrients with its 
protoplasmic network major tubes of the network become positioned in the 
non-illuminated areas (Fig. 10). 
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(a) t = h 



(b) t = 9 h 




(c) t = h 



(d) t = 9 h 



Fig. 8. Plasmodium scouting in presence of illumination obstacles, (a)-(b) experi- 
mental, (c)-(d) schemes of propagation. Oat flakes are shown by stars. 




(a) h 



(b) h 



Fig. 9. Plasmodium scouting in presence of illumination obstacles, (a) snapshot of 
experiments, (b) scheme of propagation. Oat flakes are shown by stars. 

7 Simulating Plasmodium in two-variable Oregonator 



Striking similarity in behaviour of the plasmodium and wave-fragments in 
sub-excitable media has been found in our previous papers [515] . Oregonator 
based model of the plasmodium waves impeccably matches our experimental 



results. An example of simulation is provided in Fig. 11 



At the beginning of experiments a piece of plasmodium is simulated by exciting 
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(a) 

Fig. 10. Protoplasmic network shaped by light obstacles. 
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(a) t = 200 (b) t = 400 



(c) t = 600 (d) t = 1000 



^ — I 



(e) t = 1200 (f) t = 1400 (g) t = 1600 (h) t = 1800 

Fig. 11. Simulation of Plasmodium wave diverted by recnangular illuminated do- 
main. Red and blue components of each pixel's color were defined as follows. If 
u > 0.1 then red value 255 • u, if v > 0.1 then blue value 600 • v, otherwise the 
background white color used. Time shows steps of numerical integration. 

the medium with ellipse-shaped stimulus (Fig. [Tl^t): the medium is perturbed 
by an initial excitation, when 15 x 15 sites are assigned u = 1.0 each. The per- 
turbation generates an ellipse-shaped excitation wave. Due to sub-excitability 
of the simulated medium, the initially circular wave front (Fig. [TTjp) breaks into 



two wave fragments. One wave travels North, another wave South (Fig. 11 



e). We discard South travelling wave. When the Eastern part of the North 
travelling wave enters illumination area it becomes inhibited and gradually 
disappears (Fig. [TTpf). The remaining part of the wave continues its travel as 
an individual wave, propagating North-North- West (Fig. 11 5-1). 



Simulation results for wave splitting are shown in Fig. 12 An ellipse-shaped 



zone of initial excitation (Fig. 12 1) produces two wave-fragments travelling 
North and South (Fig. fl2bc). We are concerneted with the fragment propa- 
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(a) t = 100 (b) t = 400 (c) t = 600 



A A 



(d) t = 800 (e) t = 1000 (f) t = 1200 



(g) t = 1400 (h) t = 1600 (i) t = 1800 

Fig. 12. Simulation of plasmodium wave splitting by illuminated triangle. Red and 
blue components of each pixel's color were denned as follows. If u > 0.1 then red 
value 255 • u, if v > 0.1 then blue value 600 • v, otherwise the background white color 
used. Time shows steps of numerical integration. 



gating North (Fig. 12 le). When the wave-fragment reaches illuminated tri- 
angle the part of the wave inside the triangle extinguishes due to excitability 
inhibited by light (Fig. [jjjg). In the result of the inhibition two separate wave- 
fragments are formed. One wave-fragment travels South- West another travels 
South-East (Fig.[l2}ri). 



8 Routing signals in Physarum machine 



Kolmogorov-Uspensky Machine (KUM) is denned on a labeled undirected 
graph (storage structure) with bounded degrees of nodes and bounded num- 
ber of labels [TTJMI- KUM executes the following operations on its storage 
structure: select an active node in the storage graph; specify local active zone, 
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(a) Erase(A) 



(b) Right(A) 




(c) Left(A) 



(d) Multiply(^) 



Fig. 13. Scheme of routing operations in Physarym machine. Directions of active 
zone movement are shown by arrows. Illuminated domains are rectangular and 
triangular shapes. 



i.e. the node's neighborhood; modify the active zone by adding a new node 
with the pair of edges, connecting the new node with the active node; delete 
a node with a pair of incident edges; add/delete the edge between the nodes. 
A program for KUM specifies how to replace the neighborhood of an active 
node (i.e. occupied by an active zone) with a new neighborhood, depending 
on the labels of edges connected to the active node and the labels of the nodes 
in proximity of the active node [TO] . 



Physarum machine is a biological implementation of KUM, where a node of 
the storage structure is represented by a source of nutrients (e.g. an oat flake); 
an edge connecting two nodes is a protoplasmic tube linking two sources of 
nutrients corresponding to the nodes; and, an active zone is domain of space 
(which may include food sources) occupied by a propagating pseudopodium. 
A migrating plasmodium also represents an active zone. 

In Physarum machine the computation is implemented by active zone, or sev- 
eral active zones. To make the computation process programmable one needs 
to find ways of sensible and purposeful manipulation with the active zones. 
In [7] we experimentally demonstrated how active zones can be manipulated 
by dynamical addition of nutrients. Programming with nutrients is not really 
efficient, because once the source of nutrients placed in the computing space, 
it irreversibly changes configuration of attracting fields. Light inputs allow us 
for on-line reconfiguration of obstacles and thus provide more opportunities 
for embedding complex programs in Physarum machines. 



Basic operations of routing active zones are shown in Fig. 13 
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Operation Erase(A) removes, or permanently deactivates, active zone A 



(Fig. 13 1). The operation is implemented by placing large enough - so there 
is no chance of the plasmodium wave to escape the shaded area - domain of 
illumination in front of the traveling zone A. The plasmodium wave A then 
gradually disappears. The operation is experimentally exemplified in Fig |HJ 
where few active zones are canceled by rectangles of illumination. 

Operation Left(A) rotates velocity vector of zone A on angle a, —90 < a < 0; 
for Rigth(A) < a < 90 (Fig, [ilfoc). Exact value of rotation angle a depends 
on many factors, including size of illumination domain, size of traveling zone 
A, humidity, overall illumination, and fitness of the plasmodium. Based on few 
successful experiments on rotation of plasmodium waves, we could say that 
a = 45° ± 15°. 



Operation Multiply(t4) = {Al>Ar} (Fig. 13 dc) splits traveling plasmodium 
wave A into zones Al and Ar deviated slightly left and right, comparing to 
original direction of A's travel. 



9 Discussion 



In laboratory experiments with plasmodium of Physarum polycephalum and 
in computer simulations we discovered that propagating plasmodium interacts 
with illuminated domains similarly to travelling wave-fragments in excitable 
media. We demonstrated that plasmodium waves can be diverted, annihilated 
and split by properly arranged shapes of illumination. The light-induced diver- 
sion of the plasmodium waves can be used as operations of signal and process 
routing in plasmodium-based implementations of general purpose storage ma- 
chines [1] . Future studies will concern with programming Physarum machines 
by combination of chemoattractants' gradients and illumination domains. 

Usually quite high intensity of illumination is used in studies of Physarum 
response [TT|23] . Thus in [TTJ reported threshold intensity for blue light is 
about 1,500 Lux. In our experiments the plasmodium was controlled by much 
less illuminated shapes, with maximum intensity of illumination of 50-70 Lux. 
This is why we did not observe strong photophobic reaction but rather gen- 
tle and often subtle changes in Physarum behavior. The moderate intensity 
of illumination used in our experiments shows that we can achieve purpose- 
ful behavior of the plasmodium with light stimulation without causing any 
adverse reactions. 

Geometrically shaped domains of illumination, utilised in our experiments, 
provided sharp boundaries between illuminated and non-illuminated areas. 
They did not form smooth gradients of illumination. Also size of the illumi- 
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nation domains was the same order as a 'wavelength' of propagating Plas- 
modium waves. These guaranteed that the illuminated shapes acted rather as 
reflectors, or mirrors, than sources of repelling fields. Therefore we envisage 
findings presented in the paper can be used in experimental implementation 
of collision-based computing schemes [I], where plasmodium waves represent 
quanta of information and illumination domains are used to route momentary 
wires (trajectories) of the information quanta. 
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